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Abstract 

By altering the polarization, phase and amplitude at the exit pupil, the inten- 
sity distribution near the focal plane of a Solid Immersion Lens(SIL) system 
can be changed. We have studied how the resolution and focal depth changes 
for a few particular cases. It was seen that by impinging radial polarization 
on a SIL system, we may obtain a rotational symmetric z-component of the 
focused wavefront with spot size similar to that predicted by scalar theory. 
We also observed that it was possible to manipulate the contributions from 
the homogeneous and inhomogeneous waves behind the SIL by changing the 
amplitude and phase distribution at the aperture. In this way it may be 
possible to improve both the resolution and focal depth of the system. 
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I. INTRODUCTION 



The maximum resolution achievable with conventional optical techniques is determined 
by the classical diffraction limit. The minimum optical spot diameter can be expressed as 
X/2NA, where A is the wavelength in air, and NA = nsina is the numerical aperture(n 
is the refractive index and a is the convergence semiangle). Fortunately, the diffraction 
limit can be circumvented by use of scanning near-field optical systems, where resolutions 
less than 50nm can be achieved. Unfortunately, near-field techniques have been troubled 
by low transmission efficiencies, and therefore poor signal to noise ratios. Although recent 
research have improved the transmission efficiency considerably, it is still only 1% at lOOnm 
spot size Q. Another way to increase the resolution is by application of a Solid Immersion 
Lens(SIL) [^,[3] . This method is currently not capable of the same resolution as the near- 
field techniques, but the light transmission efficiency is considerably better. The aim of the 
present paper is to describe how the polarization, phase and amplitude at the exit pupil 
influences the intensity distribution near the focal plane of a SIL. Although this is not a 
new issue in optics(see ref. and references therein), we believe it is of interest to gain 
further understanding of SIL-systems as they may become an integral part of the future 
data storage systems. Since our focusing system has high NA and the focal point is placed 
many wavelengths away from the aperture, the diffracted field near the focal plane can be 
calculated using the Debye approximation ||||. Let us consider focusing through a dielectric 
interface(see fig. [IJ). The electric field inside medium i(= 1,2) can be written as ||||: 



Where ki = 2irni/\ is the wavenumber, Sj = (s ix , s iy , s iz ) is the unit vector along a typical 
ray, fij is the solid angle formed by all the geometrical rays, T(sj) is the vector pupil 
distribution, which accounts for the polarization, phase and amplitude distribution at the 
exit pupil. We can find the electric field near the focal plane by matching the fields in the 
first and second medium at the interface, z = —d. The resulting electric field in the second 
medium becomes [Q: 

f f T(s\ s ) 

E 2 = C / lx ' ly exp[id(k 2 s 2z - k 1 s lz )}exp[ik 2 s 2z z]exp[ik 1 (s lx x + Si y y)}ds lx ds ly (2) 

J Jfli S\ z 

C is a complex constant (which will be ignored in the rest of this paper). The unit wave- 
vector is defined in spherical coordinates: 




(1) 



[sin(9i)cos((f>), sin(9i)sin((f)), cos{6i)} 



(3) 



The position vector can be written as(see also fig. [T]): 



r c = r c [sin(8 c )cos(4> c ), sin(9 c )sin((p c ), cos{6 c )} 



(4) 



This gives the following diffraction integral: 
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T{6 u 4>) = P(B u <t>)A(e u <t>) 



(6) 



where a is the convergence semiangle, P(9i,<fr) is the polarization and A(9i,4>) represents 
the amplitude and phase distribution at the exit pupil. We will for the rest of this paper 
assume that our optical system obeys the sine condition, A(9i) oc <Jcos(9i), see e.g. ref. 
II • 

k = n 2 cos(9 2 )cos(9 c ) + nisin(9i)sin(9 c )cos((p — <j> c ) (7) 

and 

* = d[n 2 cos(9 2 ) - n 1 cos(9 1 )} (8) 

represents the aberration function introduced due to the mismatch in refractive index. A 
detailed derivation of this integral was first presented in reference ||. We have presented it 
in a slightly different form to enlight the further discussion in this paper. We will consider 
two particular cases: l)Focusing of electromagnetic waves in a homogeneous media(air). In 
this case n\ — n% — 1. 2)Focusing of electromagnetic waves with a SIL, which means that 
d=0(note that n\ > n 2 ). In both cases the aberration function \1/ is identically zero 0. 



II. THE INFLUENCE OF POLARIZATION 



The state of the polarization incident on the focusing system will influence the resolution 
near the focal plane(see e.g. To discuss this question quantitatively, we must find a 

general expression for the polarization vector. We assume a incident polarization which may 
in general depend on the polar and azimuthal angle: 



b{9 1 ,(j>) 




The polarization vector can be written as |J: 

P(0i,0) = R-^LW^IL-WCRPo 



(9) 



R 



cos ((f)) sin((p) 
—sin(4>) cos(4>) 
1 



which describes the rotation of the co-ordinate system around the optical axis; 

C = 

which describes the change of polarization on propagation through the lens; 



cos(9i) sin(9\) 

1 
-sin(9\) cos{9i) 
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L Ci) 



cos(9j) —sin(9j) 

1 
sin{6j) cos{6j) 



which describes a rotation of the co-ordinate system into s and p-polarized vectors; 

t p 
i s 
t p 

which represents the transmission(Fresnel coefficients) through the plane dielectric interface. 
This results in: 



P(0i, 



a [t p cos(9 2 ) cos 2 (0) + t s sin 2 ((f))] + b[t p cos(9 2 )sin((f))cos((j)) — t s sin((f))cos((f))} 
a[t p cos(9 2 )cos((f))sin((f)) — t s sin(<p)cos(4>)} + b[t p cos(9 2 )sin 2 (4>) + t s cos 2 (<f))} 
—at p sin(9 2 )cos((j)) — bt p sin(9 2 )sin((f)) 



In the special case of x-polarized incident light (a = 1 and b=0) this matrix reduces to 

l u „„„(a \ i -L \ , 1 , 



P(0i, 



i {t p cos{9 2 ) + t a ) + -^(t p cos{9 2 ) - t s )cos(2(j)) 
\{tpCos{9 2 ) - t s )sin{2(j)) 



—t p sin(92)cos((J)) 



It is seen that different angles have different amplitudes, and that the polarization is depen- 
dent on 0. Another possibility was explored by Quabis et. al. 0, who found that radial 
polarized light may increase the resolution. We extend their analysis to a SIL system. We 
may write a(<p) = coscf) and b(<ft) = sincj): 



P(0i, 



t p cos(92)cos((p) 
tpCos(92)sin(4>) 
-t p sin(9 2 ) 



When inserted in eq. (|), this gives: 

E 2x = ill^cosifa) 



(10) 



E 2y = U[ ad sin 



(11) 



E 



2z 



rrad 



(12) 



jrad = / B{9 1 )tpSin{9 2 )sin{9 l )J {k l r c sin{9 l )sin{9 c ))exp{ik^ (13) 
Jo 



ret 

rrad _ / B{9 1 )t p cos{9 2 )sin{9 1 )J 1 {kir c sin{9 1 )sin{9 c ))exp{ik Q ^ (14) 

J 



where B{9\) = Jcos(9i), and J n is the Bessel function of the first kind, of order n. The 
z-component is completely independent of the azimuthal angle, and its importance increases 
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with increasing NA. Let us assume that we are able to find a medium which is only sensitive 
to the z-component of the polarization. That is, we may use the z-component, and disregard 
the x and y-components. Fig. |2] shows the energy density at the bottom surface(the focal 
plane) of a SIL with n% = 2(and n 2 = 1), for both radial(solid line) and x-polarized incident 
light (the dashed line represents <ft c = 0°, whereas the dotted line represents C = 90°). In this 
example, we assume that a = 60° and the wavelength Ao = 635nm. Fig. [| indicates that it 
is possible to increase the resolution, and perhaps reproduce the result predicted by scalar 
theory. However, it is worth noting that in practice one must combine the radial polarizer 
in front of the focusing system with an annular aperture, both to increase resolution and to 
avoid a singularity on the axis. This was discussed by Quabis et. al. ||. Therefore, if one 
could find a material which is sensitive only to the z-component, a more detailed calculation 
using an annular aperture is required to determine the real resolution. 



III. THE INFLUENCE OF AMPLITUDE AND PHASE 



The spatial amplitude and phase distribution(at the exit pupil) will also influence the 
energy distribution near the focal plane(see e.g. 0,0)0) • As with the incident polarization, 
there are countless ways to alter the phase and transmittance. We will limit our discussion 
to the special class of amplitude/phase filters called annular apertures(see ref. [p~0| , pTT|] and 
references therein). Annular apertures that modifies the transmittance and/or phase dis- 
tribution at the exit pupil may improve the resolution at the expense of higher sidelobes. 
In this section we will discuss the effect of such an aperture in front of a focusing system 
containing a SIL. The transmittance through an annular aperture can be expressed as: 



f Ti if < #i < a t 
T 2 if cti < Q\ < a 2 

Ti if < 6 1 < OLi 



T/v if cxn~i < 0i < cxn 

where T are complex constants(phase and amplitude) for the various zones in the aperture, 



and B(6i) = ycos(6i). If the incident light is polarized in the x-direction(a=l, b=0), the 
resulting field can be found by inserting the transmittance into eq. (pj): 



E 2x =i[Io + I 2 cos(2(j) c )} 



(15) 



E 2y = il 2 sin(2(f) c 



(16) 



E 2z = 2I 1 cos(<f> c ) 



(17) 



where 



N rcti 

I = y (Ti - T i+1 ) / 1 A^iOA Jo(A;ir c sm(0i)sm(e c ))ea;p(iA;o*)ejrp(iA;2^cos(0 2 ))^i (18) 
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N ran 

h = Y,(Ti-T i+1 ) / A lx (9 1 )J 1 {k 1 r c sin{9 1 )sin{9 c ))exp{ik ^)exp{ik 2 zcos{9 2 ))d9 1 (19) 



=1 



iV 

^(T, - T i+1 ) / A 2 %9 1 )J 2 (k 1 r c sin(9 1 )sin(9 c j)exp{ik ^)exp{ik 2 zcos{9 2 ))d9 1 (20) 



A te = B{9 x ){t s + tpoos(0 2 ))sin(0i) (21) 
A 1 * = B(9 1 )t p sin(9 2 )sin(9 l ) (22) 

A 2x = B(9x)(t a - t p cos(9 2 ))sin(9 1 ) (23) 

Where T N+1 is defined to be zero, and we have used that i+1 = Jo 1+1 — Jq* to derive these 
equations. Similar expressions can be applied to the case of radial polarized light. From 
eqs. (|IBD, ( |KSD and (PC|) we can see that the total field is a sum of electric fields from zones 
with increasing angular extent, and we expect that the terms with the smallest must be 
carefully balanced against each other in order to reduce the energy in the sidelobes. For 
high NA systems, the z-component becomes particularly important, and may increase the 
sidelobe intensity substantially. Thus it is necessary to keep the center peak ratio high in 
order to avoid a dominating z-component, which may increase the the sidelobes as well as 
the spot size. In total one must not only balance the terms from the various zones, buth 
also keep the sidelobes due to the z-(and y) component at an acceptable level. This is a 



difficult task, and is best treated by numerical analysis, e.g. by binary search methods |TT 
To keep the physics simple, we will limit ourselves to apertures with two and three zones. 
The profiles(the time-averaged energy density distributions) will be normalized, for simple 
comparison between systems. 



A. Focusing in a homogenous media 

In order to compare with the situations occuring in a SIL-system, we first observe what 
happens when we place a three-zone aperture in front of a system focusing in a homogeneous 
media. As pointed out by Ando [|HJ, the center-peak intensity ratio can be maximized by 
using a phase aperture, and the sidelobes can be made small by maintaining the same phase 
for the light passing through the center and outer portion of the aperture. On the basis of 
these results, let us assume that T\ — 1, T 2 — — 1 and T 3 = 1. Such an aperture can be 
produced by e.g. lithographical methods, and was experimentally tested by Ando et. al. 
|12]| . Let «i = 15°, a 2 = 30°, a 3 = 60° and A = 635nm. Fig. shows the focused beam 



profile for the three- zone phase aperture when C = 0° (solid line) and (j) c = 90° (dash-dotted 
line). For comparison we have also plotted the profile of the focused beam with no aperture 
when C = 0° (dashed line) and C = 90° (dotted line). The resolution due to the three-zone 
aperture is improved as compared to no aperture, but the peak sidelobe intensity is almost 
20%. We could probably bring this number down by increasing the number of zones in 
the annular aperture. Fig. [| shows the axial energy density distribution for the three-zone 
annular aperture (solid line) and no annular aperture (dashed line). Note that the focal depth 
is larger with the three-zone aperture. 
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B. Focusing with a SIL 



The last ten years much research has been done on SIL systems PJ3|,P^T5| . As mentioned 
in the introduction, a major goal has been to improve the resolution and focal depth. To 
that end, an interesting possibility is the application of an annular aperture in front of a SIL. 
As have been pointed out by Milster et al ||, plane waves incident on the bottom surface 
of the SIL will experience total reflection above the critical angle, 8 C = arcsin(ri2/ni). 
Thus we may divide the plane waves at the exit pupil into two parts; a homogeneous part 
where the angles are smaller than the critical angle at the interface, and an inhomogeneous 
part where the plane waves experience total reflection at the interface, corresponding to 
evanescent waves below the bottom surface of the SIL. Let us first see what happens when 
we put a simple two- zone annular aperture with 7\ = and T2 = 1 in front of the SIL. The 
dotted lines in figs. |5| and show the transverse profiles behind the bottom surface of the 
SIL(r2x = 2.0, n 2 = 1.0 and A = 635nm) when (fi c = 0° and C = 90°, respectively. We 
assume that oti = 54° and a 2 = 60°. For comparison, the profiles for a SIL with no annular 
aperture are also shown(the dashed lines in figs. |5| and |6|). Note that as the center (dark) disk 
increases, the profile becomes more compressed, but the sidelobes increases. The increase 
in resolution is more pronounced when C = 90° than C = 0°, since the z-component of the 
electric field depends on cos(p c (which is zero when C = 90°, see eq. (0)). Fig. |7] shows 
the axial distribution behind the SIL with(dotted line) or without (dashed line) a two-zone 
annular aperture. When the center disk increases, the focal depth decreases, contrary to 
what happens during focusing in a homogeneous media |J. This behaviour occur since 
the evanescent waves becomes more dominating when we block out the homogeneous waves, 
and since the amplitude of evanescent waves are decreasing exponentially behind the SIL. To 
increase the light transmission efficiency, one may replace the two-zone amplitude aperture 
with an axicon, thus creating a bright ring. 

Next we place a three-zone phase aperture (with the same properties as in the previous 
section) in front of the SIL. In particular, we are interested in observing what happens 
when l)«i < a 2 < 8 C and 2)a 2 > a± > 9 C . The first case(«i < a 2 < C ) is in some ways 
similar to focusing in a homogeneous media. That is, we expect the spot size to decrease, if 
the convergence semiangles have proper values. However, we also expect the focal depth to 
decrease, since we balance terms with homogeneous waves against each other, thus increasing 
the importance of the evanescent waves. As an example, let us consider the same SIL-system 
as above, but now replace the two-zone aperture with a three-zone annular aperture with 
cti = 20°, a 2 = 30° and a 3 = 60° . The dash-dotted lines in figs. ^| and |6] shows the resulting 
profiles when <fr c = 0° and <p c = 90°, respectively. Note that the profile is slightly smaller 
than the dashed line (no annular aperture) when C = 90°, whereas for <p c = 0° it is slightly 
larger. The dash-dotted line in fig. [j] shows the axial distribution. As expected, the focal 
depth is smaller with the annular aperture. 

In the second case(a 2 > cti > 9 C ) the opposite behaviour may take place. That is, we are 
able to balance terms with inhomogeneous waves against each other, thus increasing the im- 
portance of the homogeneous waves. In this way we can produce some kind of enhancement 
of the focal depth, and if the the angles a\ and a 2 are properly chosen, the beam profile may 
become narrower as well. As an example we consider an annular aperture with ot\ = 31° 
and a 2 = 43°. The solid lines in figs. || and ^ are representing the transverse profiles when 
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4> c = 0° and 0c = 90°, respectively. Now the situation is almost opposite compared with the 
previous annular aperture. That is, the profile is smaller than the dashed line(no aperture) 
when C = 0°, but larger when C = 90°. The axial distribution(solid line in fig. |?|) is 
decaying more slowly than with no annular aperture, which confirms the discussion above. 

IV. CONCLUSION 

The choice of incident polarization, phase and amplitude may change the intensity dis- 
tribution near the focal plane of a focusing system. It was seen that by impinging radial 
polarization on a SIL system, we may obtain a rotational symmetric z-component of the fo- 
cused wavefront with spot size similar to that predicted by scalar theory. We also discussed 
how annular apertures may change the resolution and focal depth. It was shown that when 
we place a three-zone phase aperture in front of a focusing system, both the resolution and 
focal depth may increase. When applied to SIL-systems, we observed that it was possible to 
manipulate the contributions from the homogeneous and inhomogeneous waves behind the 
SIL by changing the phase and transmittance distribution at the aperture. In this way it 
may be possible to improve both the resolution and focal depth. Future investigations may 
be directed towards the use of more than three zones to improve the resolution and focal 
depth. We would like to point out that only a few special cases were discussed in this paper, 
and that our examples were not optimized for any particular applications. 
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FIGURES 



FIG. 1. The general geometry for focusing through a planar dielectric interface, located at 
z=-d. The focal plane is located at z=0. 

FIG. 2. The transversal time-averaged electric energy density when n\ = 2, n 2 = 1 and a = 60°. 
l)z-component of the radial polarized focused beam. 2)Total energy density for incident x-polarized 
light(a=l and b=0) when (j) c = 0°. 3)Total energy density for incident x-polarized light(a=l and 
b=0) when 4> c = 90°. 

FIG. 3. The transversal time-averaged electric energy density with and without a three-zone 
annular aperture(ai = 15° and a 2 = 30°) when n\ = n 2 = 1 and a = 60°. l)Three-zone annular 
aperture when (j) c = 0° 2)No annular aperture when (j) c = 0°. 3)No annular aperture when <j) c = 90°. 
4) Three- zone annular aperture when <j> c = 90° 

FIG. 4. The axial time-averaged electric energy density when n± = n<i = 1 and a = = 60°. 
l)Three-zone annular aperture with a.\ = 15° and «2 = 30°. 2)No annular aperture. 

FIG. 5. The transversal time-averaged electric energy density when <p c = 0°, n\ = 2, ri2 = 1 
and a = Q3 = 60°. l)Three-zone phase aperture with ol\ = 31° and 02 = 43°. 2)No annular 
aperture. 3)Two-zone annular aperture with a\ = 54°. 4)Three-zone aperture with a\ = 20° and 
a 2 = 30°. 

FIG. 6. The transversal time-averaged electric energy density when <p c = 90°, n\ = 2, n 2 = 1 
and a = Q3 = 60°. l)Three-zone phase aperture with a.\ = 31° and a 2 = 43°. 2)No annular 
aperture. 3)Two-zone annular aperture with a\ = 54°. 4)Three-zone aperture with a\ = 20° and 
a 2 = 30°. 

FIG. 7. The axial time-averaged electric energy density when n\ = 2, n 2 = 1 and a = 03 = 60°. 
l)Three-zone phase aperture with a.\ = 31° and a 2 = 43°. 2)No annular aperture. 3)Two-zone 
annular aperture with ol\ = 54°. 4)Three-zone aperture with a.\ = 20° and a 2 = 30°. 
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Figure [| 
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